Phase structure of the two-fluid proton-neutron system 
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The phase structure of a two-fluid bosonic system is investigated. The proton-neutron interacting 
boson model (IBM-2) posesses a rich phase structure involving three control parameters and multiple 
order parameters. The surfaces of quantum phase transition between spherical, axially-symmetric 
deformed, and SU£„(3) triaxial phases are determined. 
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The phase structure of quantum many-body systems 
has in recent years been a subject of great experimental 
and theoretical interest. Models based upon algebraic 
Hamiltonians have found extensive application to the 
spectroscopy of many-body systems, including nuclei [lj 
and molecules . Applications to hadrons have also been 
developed For certain specific forms of its Hamilto- 
nian, an algebraic model exhibits dynamical symmetries. 
In the classical limit, these dynamical symmetries corre- 
spond to qualitatively distinct ground-state equilibrium 
configurations, which constitute the phases of the sys- 
tem 0, ■ The phase structure has been studied in de- 
tail for algebraic models describing systems composed of 
one species of constituent particle ("one- fluid" systems), 
in particular the interacting boson model (IBM) [l| for 
nuclei. While one-fluid systems are described by a sin- 
gle elementary Lie algebra, usually U(n), multi-fluid sys- 
tems are described by a coupling of such Lie algebras, 
Ui(n)®U a (n)®.-- 0,12. This more involved algebraic 
structure naturally leads to a richer phase structure. 

In the present work, the phase structure of a system 
comprised of two interacting fluids is investigated. The 
proton- neutron interacting boson model (IBM-2) 00, 
in which proton pairs and neutron pairs are treated as 
distinct constituents, is considered. While the one-fluid 
IBM exhibits three dynamical symmetries^separated by 
first and second order phase transitions laLDl, the IBM-2 
supports four dynamical symmetries 0,13, and thus in- 
herently has a higher-dimensional phase diagram. More- 
over, the phase structure is found to posess qualitatively 
new features. Due to the complexity of the problem, 
a combination of analytic and numerical methods have 
been applied in this analysis. Preliminary investigations 
of the IBM-2 phase structure have been presented in 
Refs. H1E3- 

Before proceeding, let us briefly summarize the IBM-2 
Hamiltonian and the dynamical symmetries it supports. 
Operators in the IBM-2 are constructed from the gen- 
erators of the group U w (6) ® U„(6), realized in terms of 
the boson creation operators sj,,o and d) p p (where p rep- 
resents 7T or v, and fi——2, . . . , 2) and their associated 



annihilation operators, acting on a basis of good boson 
numbers N n and N v . A schematic Hamiltonian which 
retains the essential features of the model is the F-spin 
invariant Hamiltonian (e.g., Rcf. 13]) 

H = e(h dTr + n dv ) + + Q*") • + QX^ (r) 



where n d p = d/-d p and Q 
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d p )( 2 \ It is convenient to also introduce "scalar" and 
"vector" parameters Xs=\{x^+Xv) and Xv=\{x-k~Xu)- 
Three of the IBM-2 symmetries occur for Yy=0 and 
have direct analogues in the one-fluid IBM [8j: U wv (5) 
(k = Q) , for which the geometric interpretation is that of 
undeformed proton and neutron fluids, SO,n/(6) (e = 0, 
Xir — Xv = 0), yielding deformed, 7- unstable structure, 
and SU,r^(3) (e = 0, x-* — Xv = — s/7/2), for which 
prolate axially symmetric structure is obtained. [The 
complementary case X7r = Xr/=+\/7/2, giving oblate axi- 
ally symmetric structure, is distinguished by the notation 
SU,r„(3).] However, a symmetry special to the IBM-2, 
denoted SU*„(3), is obtained for e — 0, x-k — ~V7 /2, and 
Xv = + V7/2 The equilibrium configuration consists 
of a proton fluid with axially symmetric prolate deforma- 
tion coupled to a neutron fluid with axially symmetric 
oblate deformation, with their symmetry axes orthogo- 
nal to each other [flL llilllij . yielding an overall composite 
nuclear shape with triaxial deformation. To avoid ambi- 
guity, we shall adopt here the notation 811*^(3) for the 
complementary case x?r = +V7/2 and Xv = —y/7/2, for 
which the proton and neutron deformations are inter- 
changed. 

The classical limit of the IBM-2 is obtained by eval- 
uation of the expectation value of H for the coherent 
state K, at 2) ; N v , ai 2) )cx (4 >0 + £ M aj&dt^)*- (4, + 
J2/j, a< v}idl, l j,) Nl '\0)- This yields an energy surface £=(H) 

as a function of the coherent state parameters a p ,p,. 
(2) 

The Op/fj, are interpreted geometrically as the quadrupole 
shape variables [lfj for the proton and neutron fluids 
and are equivalent to four deformation parameters (/3 n , 
j n , (3 U , and j u ) and six Euler angles (d\ v , 8 27r , 03n, Q\v, 
02v, and #3„). By rotational invariance, £ depends only 
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upon the relative Euler angles between the proton 
and neutron fluid intrinsic systems, not the 9^ and 9 lv 
separately. Minimization of £ {(3^, 7^, V , j y , #1, i?2j $3) 
with respect to its parameters yields the classical equilib- 
rium configuration of the ground state. The terms (fid-jr), 
(ndv), (Qtt ■ Qtt), and (Q v ■ Q w ) contributing to £ involve 



only a single fluid and arc thus known from the IBM (sec 
Ref. The expectation value (Q^ ■ Q v ) of the inter- 

action term is obtained by the methods of Refs. 0, 
as a function of all seven possible parameters ((3^, 7^ , /?„, 
Ju, #1, i?2i and i? 3 ), 



■ = TTTW^TTW^ * + + ^ * x ^ 2) ) (2) 1 ' ^ * + ^ + *>lf® * x ( 2 ) 



where 

afl = P P cos lp Dl*(9 lp , 9 2p , 9 3p ) 
+ -Lp p sm lp [Dl;(9 lp , 9 2p , e 3p ) + Dl*_ 2 (0 lp , 9 2p , 9 3p )] 

(3) 

and where the Euler angles may be chosen to be 9^ = 
and 9i v = i?j by rotational invariance. The (hd P ) are 
linear in or N v , while the (Q p ■ Q p >) are quadratic. It 
is thus convenient to reparametrize the Hamiltonian 

as 

H = + n dv ) - Jl(<&- + Ql" ) ■ (Qx- + 

(4) 

where N = N„ + N v , so that the energy function £ is 
independent of N at fixed ratio N^/N^. This definition 
also condenses the full range of possible ratios e/n onto 
the finite interval < £' < 1. An overall normalization 
parameter for H has been discarded as irrelevant to the 
structure of the energy surface. There are thus three con- 
trol parameters — x-m and \ v — for this Hamiltonian. 

A simple categorization of the possible Euler angle and 
7 P values for the equilibrium configurations for certain 
IBM- 2 energy surfaces has been presented in Ref. [l5j . 
For a class of Hamiltonians including the present one |@J, 
it is found that the global minimum only occurs for van- 
ishing relative Euler angles, i.e., for aligned proton and 
neutron intrinsic frames. This effectively reduces the 
number of order parameters for the system from seven 
to four — 7tt, Pv, and 7„. 

The orders of phase transitions are, in the present 
study, determined according to the Ehrenfest classifica- 
tion: a phase transition is first order if the first deriva- 
tive of the system's energy is discontinuous with respect 
to the control parameter being varied, second order if 
the second derivative is discontinuous, etc. Where the 
system's energy is obtained, as in the present classical 
analysis, as the global minimum of an energy function £ , 
a first order transition is usually associated with a discon- 
tinuous jump in the equilibrium coordinates ("order pa- 
rameters") between distinct competing minima. Second 



or higher order transitions are associated instead with a 
continuous evolution of the equilibrium coordinates, as 
when an initially solitary global minimum becomes un- 
stable (posessing a vanishing second derivative with re- 
spect to some coordinate) and evolves into two or more 
minima. It should be noted that, whenever the order 
of a phase transition is obtained by numerical analysis, 
application of the Ehrenfest criterion is limited by the 
ability to numerically resolve sufficiently small disconti- 
nuities, especially a consideration for points of first-order 
transition very close to a point of second-order transi- 
tion. Moreover, problems with the classification scheme, 
not addressed here, may arise at the boundaries of the 
parameter space or when the Hamiltonian posesses addi- 
tional symmetries. 




FIG. 1: Phase diagram for the SU™(3)-SO,r„(6)-SU;„(3) 
plane (£' = 1) in the parameter space of the Hamiltonian 
of |@J, showing the second-order phase transition curves © 
for the cases N^/N^ = 1 (solid) and N^/N v = 4 (dashed). 
The diagram is rotated to allow more direct comparison with 
Fig.H 
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We begin our analysis with an analytic study of the 
phase structure of the Hamiltonian for £' = 1, which 
encompasses the SO WI/ (6), SU 7ri/ (3), and SU*„(3) dynam- 
ical symmetries (see Fig. The global equilibrium in 
this case is always deformed (f3„ > and (3 V > 0). Sur- 
rounding the SU7rj/(3) dynamical symmetry is a region 
of parameter space in which the equilibrium deforma- 
tions are axially symmetric (7^ = 7 „ = 0°), and a sim- 
ilar region surrounds the SU,n,(3) dynamical symmetry 
(7tt = Jv = 60°). Taking the SU 7ri y(3)-like region for speci- 
ficity, the global minimum occurs for 
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(5) 



At the boundary of this region, axial equilibrium defor- 
mation gives way to triaxial deformation, with 7^ and/or 
7„ nonzero. This transition occurs continuously on the 
locus of points at which the minimum given by © first 
becomes unstable with respect to 7 deformation. Since 
£ depends upon both 7ir and j u , instability occurs when 
the directional second derivative of £ first vanishes along 
some "direction" in (7^, 7 „) coordinate space, which may 
generally happen before either d 2 £/d~f 2 or d 2 £/d^ 2 van- 
ishes individually. The equation describing the bound- 
ary curve in x-k and Xu is most compactly expressed in 
terms of the corresponding equilibrium values (3^ and [3 V 
from © as 

, 9&&G#-l) + &(2/£-l)(# + l) 



-2) 2 

^-1)+ ^(2^-i)o? + i) 



2/3,(^-2)2 



(6) 



This curve is shown in Fig. ^ Along the SU,r„(3)- 
SU*„(3) line (x-k = -V7/2), for N w /N u = 1 the transition 
occurs at Xv ~ 0.4035, at which point the global minimum 
becomes soft with respect to 7^ at fixed 7^ = 0°. 

Returning to the full parameter space for the three- 
parameter Hamiltonian of (gj , limited but useful analytic 
results can also be obtained for the transition between 
undeformed (f3 n = (3 V = 0) and deformed structure. The 
parameter space is illustrated in Fig. [21 

For xv=0, *-e., in the "base" plane in Fig. El the anal- 
ysis is closely related to that for the one-fluid IBM. The 
equilibrium configurations all have @ w = (3 V (= (3) and 
are identical to those obtained for the IBM Hamiltonian 
fliBM = [(1 - C)/N]h d - (i'/N 2 )& ■ With Y = X S, 
the phase structure of which is well known [Sj, A 
second-order transition between undeformed (j3 = 0) and 
deformed {(3 7^ 0) structure occurs at the parameter val- 
ues £' = 1/5 and x = 0, for which the minimum in the 
energy surface at (3 = is unstable. This point lies on a 
trajectory of first-order transition points, at which a dis- 
tinct minimum with nonzero (3 preempts that with f3 = 
as global minimum. 
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FIG. 2: Phase diagram of the proton-neutron interacting bo- 
son model (IBM-2) for the Hamiltonian of @, as obtained 
by numerical minimization of £ , for N^/N^ — 1. The sur- 
faces of first-order (red) and second-order (blue) transition 
between regions of undeformed, axially symmetric deformed, 
and triaxially deformed equilibria are shown. Only one "quad- 
rant" of parameter space (0 < £' < 1, —y/7/2 < xs < 0, and 
< Xv < V7 /2) is included in this plot, since the others may 
be obtained by reflection. The xs and xv axes are scaled by 
so as to converge to a point at the U,r„(5) limit. 



A second-order transition from undeformed to de- 
formed structure occurs when the minimum of £ at 
(3^ = [3 U = § becomes unstable with respect to (3 de- 
formation, provided this minimum is the global mini- 
mum (that is, provided it has not been rendered irrel- 
evant by a prior first-order transition to another, de- 
formed minimum) . The derivative indicating such (3 soft- 
ness is the directional second derivative d 2 £/d0 2 along 
a ray (3^ — u^(3 and /3 U — u v f3 (i.e., fixed f3„/f3 v ) at 
fixed 7-n- and 7 „, evaluated at f3 = 0. This quantity 
is independent of xs and xv [e.g., for N^/N v = 1, 
d 2 £/dp 2 1/3=0 = (1 - 3£')(t4 + u 2 u ) -AguvUy cos( 7t - 7 „)]. 
The minimum at zero deformation first becomes unsta- 
ble at £' = 1/5, where it is soft against deformations with 
Ar = f3 v and 7^ = j v (= 7), for any value of 7. 

The curve of first-order phase transition in the plane 
Xv = arises from competition between the undeformed 
minimum and one with j3 n = j3 v 7^ and 7^ = 7„ = 0°. 
The special "slice" £(/3„ = [3,(3 V = /3,j v = 0,7„ = 0) of 
the energy function, which includes both these minima, 
is found for N^/N^ — 1 to be independent of xv at fixed 
£' and xs> invariant along any vertical line in Fig. [21 
Thus, the first-order phase transitions occuring in the 
plane xv = at £' < 1/5 "propagate" out of this plane. 
[To the approximation that (3^ f=a (3 V for the deformed 
minimum, the first-order transition surface for xv 7^ is 
obtained by vertical extension of the one-fluid IBM tran- 
sition trajectory in Fig. [5J] The occurence of a second- 
order phase transition at £' = 1/5 is thus precluded ev- 
erywhere except along a vertical line in parameter space 
extending through the one-fluid IBM second-order tran- 
sition point. It is verified numerically, as described be- 
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low, that the undeformed minimum is indeed global along 
this line. The line = 1/5 and xs = is thus a locus of 
second-order transition between zero and nonzero defor- 
mation. These results are readily generalized to arbitrary 
Nn/Nv, for which the invariance off occurs along lines 
of constant N^X-n + N„Xv and the line of second-order 
phase transition obeys X-n/Xu = —N v /N v . 

The remainder of the phase diagram is obtained by nu- 
merical minimization of the energy surface with respect 
to Pit, 77r, f} v , and 7„. For robust identification of the 
global minimum, £ is first evaluated at each point on a 
fine mesh in these coordinates (A/3 p = 0.03, A7 P = 2°), 
and all points which are discrete local minima of £ rel- 
ative to the neighboring mesh points are identified. The 
0n, 7tt, Pv, and ^ v values for these minima are then re- 
fined by an iterative method. The global minimum is 
identified from among these. 

The phase diagram obtained in this fashion is shown in 
Fig. [3 for the case N^/N^ = 1. Numerical investigation 
of the behavior of £ at the transition points allows the 
Ehrcnfcst criterion to be applied, and it appears that the 
axial-triaxial transition is everywhere second order. 

The IBM-2 phase diagram obtained here provides a 
framework for studying the transition between axial and 
triaxial structure in nuclei. Triaxial nuclear deformation 
might arise from several different sources: higher-order 
(cubic, quartic, etc.) interactions in an essentially one- 
fluid nucleus [17| , distinct deformations of the proton and 
neutron fluids as discussed here, or the presence of config- 
urations involving hexadecapole nucleon pairs The 
present analysis provides insight into the conditions un- 
der which proton-neutron triaxial deformation may occur 
and the nature of the transition to such structure. Phe- 
nomenological studies extending this work should make 
use of a more realistic Hamiltonian involving a Majo- 
rana contribution and different strengths for the Q p ■ Q p i 
terms {e.g., Ref. |l3). Such studies will provide guidance 
for the experimental investigation of triaxial structure as 
progressively more neutron rich nuclei become accessible. 

The present analysis may serve as a model for the 
study of other multi- fluid bosonic systems. Within nu- 
clear physics, the U corc (6) <g> U s ki n (6) description of core- 
skin collective modes in neutron rich nuclei can be 
treated similarly. In molecular physics, a study of the 
phase structure is in progress for the vibron model with 
two vibronic species [2j, where it is relevant to coupled 
vibronic bending modes in acetylene. Another area of 
potential application of the method is to atomic Bose- 
Einstein condensates. Scissors modes, introduced orig- 
inally within the framework of nuclear two-rotor mod- 
els and the IBM-2 [2l| , have been observed for oscil- 
lations of a single-constituent Bose-Einstein condensate 
relative to an anisotropic potential |22j . Experiments are 
planned to produce condensates of two different atomic 
species and to study the scissors modes between them. 
Triaxial deformations of the type considered here could 



then occur, and a directly analogous analysis could be 
applied. 

The exotic features of the phase diagram considered 
here, arising from the presence of multiple control param- 
eters and multiple order parameters, are likely to be en- 
countered for other multi-fluid systems as well. They will 
require a classification scheme beyond the simple Ehren- 
fest or one-parameter Landau models for their proper 
description. 
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